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ABSTRACT 

We present deep u*-, and g-band images taken with the MegaCam on the 3.6 m Canada-France-Hawaii 
Telescope (CFHT) to support the extragalactic component of the Spitzer First Look Survey (hereafter, FLS). In 
this paper we outline the observations, present source catalogs and characterize the completeness, reliability, 
astrometric accuracy and number counts of this dataset. In the central 1 deg^ region of the FLS, we reach 
depths of g ^26.5 mag, and u* ^26.2 mag (AB magnitude, Scr detection over a 3" aperture) with ^4 hours of 
exposure time for each filter For the entire FLS region (^^5 deg^ coverage), we obtained M*-band images to the 
shallower depth of u* =25.0-25.4 mag (5cr, 3" aperture). The average seeing of the observations is 0."85 for the 
central field, and ^L"00 for the other fields. Astrometric calibration of the fields yields an absolute astrometric 
accuracy of 0."15 when matched with the SDSS point sources between 18< g <22. Source catalogs have been 
created using SExtractor. The catalogs are 50% complete and greater than 99.3% reliable down to § ~ 26.5 
mag and u* ~ 26.2 mag for the central 1 deg^ field. In the shallower M*-band images, the catalogs are 50% 
complete and 98.2% reliable down to 24.8-25.4 mag. These images and source catalogs will serve as a useful 
resource for studying the galaxy evolution using the FLS data. 

Subject headings: catalogs — surveys — galaxies:photometry 



1. INTRODUCTION 

Since its launch in August 2003, the Spitzer space telescope 
has opened new IR observing windows to the universe and 
probed to unprecedented depths. The first scientific obser- 
vation undertaken with the Spitzer after its in-orbit-checkout 
period was the Spitzer First Look Survey. The survey pro- 
vided the first look of the mid-infrared (MIR) sky with deeper 
sensitivities than previous systematic large surveys^, allowing 
future users to gauge Spitzer sensitivities. For the extragalac- 
tic component of the First Look Survey (FLS), the ^4.3 deg^ 
field centered at RA=17''18'"00', DEC=H-59°30'00" was ob- 
served with the InfraRed Array Camera (IRAC, 3.6, 4.5, 5.8, 
8.0/im; Fazio et al. 2004) and the Multiband Imaging Pho- 
tometer for Spitzer (MIPS, 24, 70, 160/^m; Rieke et al. 2004). 
The main survey field has flux limits of 20, 25, 100, and 
100/iJy at wavelengths of IRAC 3.6, 4.5, 5.8, and 8.0^m re- 
spectively. The deeper verification field, covering the central 
900 arcmin^ field, has sensitivities of 10, 10, 30, and 30/iJy 
at IRAC wavelengths. (Lacy et al. 2005) Also, the FLS veri- 
fication strip field has 3(7 flux limits of 90/iJy at MIPS 24/im 
(Yan et al. 2004a), and 9, and 60 mJy at 70, 160^m (Frayer et 
al. 2005). The reduced images and source catalogs have been 
released (Lacy et al. 2005; Frayer et al. 2005). 

In order to support the FLS, various ground-based ancil- 
lary datasets have also been obtained. A deep /?-band sur- 
vey to RAB{5a)=25.5 mag (KPNO 4 m, Fadda et al. 2004) 
and a radio survey at L4GHz to 90^Jy (VLA, Condon et 
al. 2003) have been performed to identify sources at differ- 
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ent wavelengths. Imaging surveys with the Palomar Large 
Format Camera (LFC), g' , r' , i' and z' band data covering 1- 
4 deg^ have also been taken since 2001. The Sloan Digital 
Sky Survey also covers the FLS field. Moreover, spectro- 
scopic redshifts have been obtained for many sources using 
the DEep Imaging Multi-Object Spectrograph (DEIMOS) at 
the Keck observatory as well as Hydra instrument on the Wis- 
consin Indiana Yale NOAO (WIYN) observatory. Utilizing 
all the Spitzer and the ground-based ancillary datasets, many 
science results have come out already from the FLS. The FLS 
studies range from the infrared source counts to the proper- 
ties of obscured galaxies such as submm galaxies, red Active 
Galactic Nuclei(AGNs), and Extremely Red Objects (e.g. Yan 
et al. 2004b; Fang et al. 2004; Marleau et al. 2004; Frayer et 
al. 2004; Appleton et al. 2004; Lacy et al. 2004; Choi et al. 
2005). 

In order to broaden the scientific scope of the FLS, we have 
obtained new ground-based ancillary datasets, deep M*-band 
and g-band imaging data that will be presented in this paper. 
The exploration of the cosmic star formation history at z < 3 
is one of the key scientific issues that the FLS is designed to 
address. While studying the infrared emission is an extremely 
valuable way to understand the obscured star formation his- 
tory of the universe, a complete census of the cosmic star for- 
mation history requires the measurement of unobscured star 
formation as well. Such a measurement can be done effec- 
tively using the rest-frame ultraviolet (UV) continuum, which 
represents either instantaneous star formation activity of mas- 
sive stars (Far UV; FUV), or traces the star formation from in- 
termediate age stars (Neai- UV; NUV) (Kennicutt 1998). The 
FUV and NUV emission of galaxies at z < 1 redshifts to the 
M-band; therefore, the addition of the M-band enables us to de- 
scribe the star formation history at z < 1 to the full extent. 
Also, using the deep M-band and g-band data, it is possible 
to select M-band dropout galaxies - star forming galaxies at 
z ~ 3. Especially of interest are bright, massive Lyman Break 
Galaxies which are rare, but important to constrain galaxy 
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TABLE 1 

Characteristics of the filters used in MegaCam observations 



Filter 


u* 


g 


central wavelength (A) 


3740 


4870 


wavelength range 


3370-4110 


4140 - 5590 


mean transmission (%) 


69.7 


84.6 



evolution models since they can tell us how massive galax- 
ies are assembled in the early universe. The wide area cov- 
erage of the FLS allows us to select such rare, bright LBGs. 
However, for the selection of LBGs, it is necessary to obtain 
M-band data and g-band data that matches the depth of the 
Spitzer data. In addition, the M-band and g-band datasets can 
be used to test extinction correction methods, such as the es- 
timation of dust extinction based on the measurement of UV 
slope (Calzetti et al. 1994; Meurer et al. 1999; Adelberger & 
Steidel, 2000). When used together with other ground-liased 
ancillary data, our u- and g-band data will be very helpful for 
improving the accuracy of photometric redshifts. 

With these scientific applications in mind, we present the 
u*- and g-band optical observations made with the MegaCam 
at Canada France Hawaii Telescope (CFHT). The dataset is 
composed of i) deep u*-, g-band data for the central 1 deg^ of 
the FLS, and ii) M*-band data for the whole FLS. In §2, we 
describe how our observations were made. The reduction and 
calibration procedures are in §3. Production of source cata- 
logs is described in §4. Finally in §5, we show the properties 
of the final images and the extracted catalogs. Throughout 
this paper, we use AB magnitudes (Oke 1974). 

2. OBSERVATIONS 

The observations were made using the MegaCam on the 3.6 
m Canada France Hawaii Telescope (CFHT). MegaCam con- 
sists of 36 2048x4612 pixel CCDs, covering 0.96 degxO.94 
deg field of view with a resolution of 0.187" per pixel (See 
the detailed description in Boulade et al. 2003). The avail- 
able broad-band filters are similar to SDSS ugriz filters, but 
not exactly the same especially in the case of M*-band. We 
used u* and g filters for our observations. The filter character- 
istics are summarized in Table 1, and their response curves 
are compared to the SDSS u- and ^-band filters in Figure 
1. Note that these response curves include the quantum effi- 
ciency of the CCDs. The response curves of the two g-bands 
are nearly identical in their shapes, but the M*-band response 
curve shows a different overall shape compared to the SDSS 
M-band, in such a sense that u* filter is redder than SDSS u. 
The difference leads to as much as a 0.6 mag difference be- 
tween the CFHT & SDSS M-magnitudes. This will be investi- 
gated further in the photometry section (§5.2). 

Our images on the FLS field were taken with the queued ob- 
servation mode in two separate runs (03B and 04A semesters). 
The first observing run occurred in August 23-29 (M*-band) 
and September 19-22 (g-band), 2003. During this 03B run, 
we observed the central 1 deg^ of the FLS field, which cov- 
ers the entire FLS verification strip. Images are taken with 
medium dither steps of about 30", at 12 (m*) or 20 (g) differ- 
ent dither positions. The exposure time at each position was 
1 160 seconds for M*-band, 680 seconds forg-band. Therefore, 
the total integration time per pixel amounts to ~3.8 hours for 
both bands. The yellow-colored box in Figure 2, shows the lo- 
cation of this deep u*- and g-band data superimposed on the 
/?-band images (Fadda et al. 2004). 
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Fig. 1. — The response curves of CFHT u", g filters mul- 
tiplied by the telescope (mirror, optics) and CCD efficiency (solid). 
The CFHT data points are taken from a CFHT-related homepage 
(http://astrowww.phys.uvic.ca/grads/gwyn/cfhtls). The SDSS u-, g-band 
transmission curves are given as dashed lines (http://www.sdss.org). The 
CFHT filters are sUghtly redder than their SDSS counterparts. 

In the second observing run during 04A period, we obtained 
M*-band images covering the whole FLS field (^5 deg^) at a 
shallower depth. Note that we did not take g-band data be- 
cause the entire FLS is covered to a shallower depth by the ex- 
isting g-band data taken with the Palomar 5 m telescope. The 
CFHT observation was at first designed to take 10 dithered 
images at 7 different locations of the FLS area, which are 
named as FLS11-FLS32 (Table 2). The cyan-colored boxes 
marked with the dashed line in Figure 2 indicate FLSll- 
FLS32 fields. The dates when the observations made span 
quite a wide range from April to July 2004. Some of the fields 
have all dithered images observed on a photometric condition 
(FLS21, FLS22 & FLS31), but other fields contain the im- 
ages that are non-photometric. Due to the bad weather, we 
lost some of the queued observing time and several observed 
images were not validated. Each dithered image was taken 
with the exposure time of 680 seconds. Although the origi- 
nal plan was to achieve 1.9 hours of total integration at each 
location, there is a field-to-field variation on the image depth 
between ^1 hour and ^2 hours. Table 2 gives the summary 
of the 03B and 04A observations. The seeing condition of the 
whole observation sessions ranged from 0"85 to 1"08. 

Between August 2003 and September 2003, one of the 36 
MegaCam CCDs was malfunctioning. As a result, g-band im- 
ages lack [CCD03] part of the observed field. Hence, the ef- 
fective survey area in coadded mosaic image is ~0.94 deg^ 
for the central field. 

3. DATA REDUCTION 

3.1. Preprocessing with Elixir 

AH MegaCam data obtained at CFHT are preprocessed 
through the EUxir pipeUne^, the standard pipeline for the ba- 
sic reduction of the MegaCam data, provided by the CFHT. 
Through the Elixir pipeline, the raw object frames are bias 
subtracted, flat fielded, and photometrically calibrated. In the 
EUxir photometric caUbration, the nightly photometric zero- 
point is estimated using the standard star data on that night. 
Basic astrometric calibrations are done using ~50-60 USNO 

* See Elixir website http://www.cfht.hawaii.edu/Instruments/Elixrr 



Fig. 2. — Area coverages of various datasets over the FLS. The small yellow square in the center represents the Spitzer FLS verification strip, and the big 
yellow square shows our CFHT observations in the central 1 deg^. The red/white squares correspond to the whole FLS area (IRAC/MIPS respectively) which 
covers ^4.3 deg^. The cyan squares indicate the area covered by our CFHT w*-band observations in 2004. All these marks are overlayed on the KPNO i?-band 
images. 



TABLE 2 
Observation SUMMARY 



Field ID 


Q (J2000) 


5 (J2000) 


Observation Dates 


Filter 


Total Exposure Time (sec) 


Depth (AB mag)" 


Seeing ("t 


Zeropoint'^ 


central 


17:17:01 


+59:45:08 


2003/08/25, 2003/08/29 


u* 


1160 X 12 


26.41 


0.85 


32.904 


central 


17:17:01 


+59:45:08 


2003/09/19, 2003/09/21, 2003/09/22 


S 


680 X 20 


26.72 


0.85 


33.430 


FLS 11 


17:22:09 


+60:14:47 


2004/04/29, 2004/07/12 


Ll* 


680 X 10 


25.83 


0.98 


32.226 


FLS 12 


17:14:44 


+60:14:50 


2004/06/18, 2004/07/17 


U* 


680 X 9 


25.85 


1.06 


32.186 


FLS21 


17:22:14 


+59:30:00 


2004/07/10, 2004/07/21 


U* 


680 X 9 


25.91 


0.98 


32.236 


FLS22 


17:14:47 


+59:44:45 


2004/07/10 


u* 


680 X 5 


25.57 


1.08 


32.213 


FLS23 


17:07:09 


+59:48:26 


2004/07/11 


u* 


680 X 5 


25.43 


0.89 


32.221 


FLS31 


17:20:46 


+58:49:09 


2004/07/11, 2004/07/22 


u* 


680 X 8 


25.68 


0.90 


32.188 


FLS32 


17:13:32 


+58:57:46 


2004/07/12 


Ll* 


680 X 7 


25.67 


0.95 


32.139 



^The depth of a field is calculated as 5a flux over a 3" aperture. 

''The seeing of the image was determined through a visual inspection on individual stars. 

'^The zeropoints are calibrated from original Elixir photometric solution (See §3.5). The values are given as zeropoints according to DN. 
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Fig. 3. — The figure shows the improvement in astrometric solutions after the astrometric calibration with SDSS sources. The vectors indicate the difference 
between reference SDSS positions vs detected positions. The size of the vectors are exaggerated by a factor of 250. The head of the arrows are pointing to SDSS 
positions. After the calibration, the size of the vectors remarkably decreased (more than a factor of 5) and the distortion is reduced significantly. The figure 
represents 6 adjacent CCD chips. The borders of the CCD chips are shown as dashed lines. 
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(United States Naval Observatory catalog; Zacharias et al. 
2000) and HST GSC (Hubble Space Telescope Guide Star 
Catalog; Morrison et al. 2001) reference stars per chip. In 
this calibration, the average pixel scale and the image rotation 
are derived and written as WCS keywords to the header of the 
raw image. Since the Elixir astrometric solution uses only a 
first order fit, the absolute astrometric accuracy of the Elixir- 
processed images is about 0"5-l" rms with respect to the ref- 
erence stars above. Final processed images are delivered to us 
as Multiple Extension Fits files that can be manipulated using 
IRAF' mscred. 

3.2. Post-Elixir Processing before Mosaicking 

The dehvered images are processed before mosaicking, us- 
ing various IRAF packages and tasks. The post-EUxir image 
processing includes: (1) identification of saturated pixels and 
bleed trails; (2) creation of new bad pixel masks; (3) removal 
of satellite trails by inspecting each image frame; (4) replace- 
ment of bad pixel values; and (5) removal of cosmic rays. 

As the first stage of the post-Elixir image processing, we 
identified saturated pixels and bleed trails using the ccdproc 
task in IRAF. The saturated pixels and bleed trails, which 
are the pixels showing nonlinear behavior, should be identi- 
fied before stacking images because they can affect surround- 
ing pixels during the projection of the image. Those pixels 
with values above the saturation level, 64000 DN for M*-band 
and 62000 DN for g-band, were identified as saturated pixels. 
Neighboring pixels within a distance of 5 pixels along lines or 
columns of the selected saturation pixels are also identified as 
saturated. We found that the pixel saturation occurs typically 
at the central part of the stars brighter than u* ~20 mag and 
g ~21 mag. Bleed trails, which result from charge spillage 
from a CCD pixel above its capacity, tend to run down the 
columns. We identify bleed trails as more than 20 connected 
pixels that are 5000 counts above the mean value along the 
whole image (i.e., setting the ccdproc threshold parameter 
as "mean-i-5000"). After the identification of saturated pix- 
els and bleed trails, this information is stored as a saturation 
mask. 

Next, we used the IRAF task imcombine to make new 
bad pixel masks by combining the saturation masks with the 
original CFHT bad pixel masks. The resultant masks have 
pixel values of for good pixels, greater than for bad pixels. 
In the bad pixel masks, we also indicate the location of satel- 
lite trails. To include satellite trails in bad pixel masks, we 
identified satellite trails through the visual inspection by flag- 
ging their start/end (x, y) image coordinates and widths on 
each CCD chip. The flagged rectangular line is then marked 
as bad pixels in the bad pixel masks. 

With the newly constructed bad pixel masks, we fixed the 
bad pixels using the IRAF f ixpix task. The value of the 
pixel marked in the bad pixel masks was replaced by the lin- 
ear interpolation value from adjacent pixels. The next step is 
to identify the pixels affected by cosmic rays. We applied a 
robust algorithm based on a variation of Laplacian edge de- 
tection with the program, la_cosmic (van Dokkum 2001) 
on each CCD frame. This algorithm identifies cosmic rays of 
arbitrary sizes and shapes, and our visual inspection of this 
procedure confirms effective removal of cosmic rays. After 

' IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 



the cosmic ray cleaning, we improved the astrometric solu- 
tion, which is described in the next section. 

3.3. Astrometric Calibration 

The Elixir data pipeline performs a rough astrometric cal- 
ibration on each image, and the existing header WCS key- 
words are updated with the improved astrometric solution. As 
mentioned in Elixir section (§3.1), Elixir calculates the astro- 
metric solution for each CCD chip. An initial guess is de- 
termined for the image coordinate based on the initial header 
WCS keywords, and Elixir uses HST Guide Star Catalogs or 
USNO database to calculate an astrometric solution. The so- 
lution derived by Elixir is a first order fit that gives the rms 
scatter of 0"5^1". Thus, the solution is not good enough for 
follow-up observations that require very accurate position in- 
formation such as the multi-object spectroscopy. The accurate 
astrometry is also important when combining and mosaicking 
images, since misalignment of images can lead to a loss of 
signal in the combined image. 

To improve the astrometric solutions, we adopted SDSS 
sources as reference points, and derived the astrometric solu- 
tion using msctpeak task in mscf inder, which is a sub- 
package of the IRAF mosaic reduction package, mscred. 
The SDSS sources brighter than r ~20 mag have an absolute 
astrometric accuracy of 0"045 rms with respect to USNO cat- 
alogs (Zacharias et al. 2000) and 0"075 rms against Tycho-2 
catalogs (Hog et al. 2000) with an additional 0f'02~0f'03 
systematic error in both cases (Pier et al. 2003). We also con- 
sidered using USNO and GSC catalogs; unfortunately, many 
of the USNO/GSC stars were saturated in our images, making 
it difficult to use them for astrometric calibration. 

There are other tasks such as msczero and msccmatch 
for deriving the astrometric solution, but we settled on us- 
ing msctpeak. Solutions from tasks such as msccmatch 
- which works only in TAN projection - are reflected in the 
value of CDi_j, CRPIXi header keywords. As a result, 
the astrometric solutions from such tasks have limitations and 
they do not improve the astrometric accuracy significantly. On 
the other hand, msctpeak creates a separate WCS database 
including a higher order fit solution from TNX sky projection 
geometry, and they are written as additional keywords such as 
WATi_00 j in the FITS header using mscsetwcs task. The 
resultant astrometry of the image improved significantly over 
the original Elixir astrometry solution. 

Ultimately, we used SDSS sources with 16 < m < 22 mag 
and 17 < g < 22 mag to calculate the astrometric solutions 
for the u- and g-band images, respectively. The astrometric 
calibration objects were inspected by eye in order to exclude 
heavily saturated sources. The sources used for the astromet- 
ric calibration are not restricted to stars, since the SDSS fil- 
ters and CFHT u* and g filters are similar enough that the 
centroiding problem is not a great concern. Including both 
stars and galaxies, we could get ~90 objects per chip to de- 
rive the astrometric solution. The solutions are calculated in- 
teractively for each chip, therefore there are 36 solutions for 
each dither image. Functions used for the fitting are polyno- 
mials of order 3~4, and the fitting errors are typically below 
0"1. How the astrometric solutions fair with other data are 
described in a separate section (§5.1). Figure 3 demonstrates 
the improvement of astrometric solutions after executing the 
above procedure. 

3.4. Production of Mosaic Images 
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Fig. 4. — The magnitude difference of matched objects from different 8 i<*-band fields before zeropoint correction is shown. Note that the objects lying in 
FLS21, FLS22, and FLS31 fields are brighter than the objects lying in the central field before zeropoint correction. The mean value of the magnitude offsets are 
(FLSll-central)~0.006, (FLS12-central)~0.018, (FLS21-central)~-0.028, (FLS22-central)~-0.029, (FLS22-FLS23)~-0.009, (FLS3 l-central)~-0.035, 
and (FLS32-central)~-0.001. Since the magnitude offsets of 1^^.821, 1'LS22, and ^..831 fields are consistent within 0.01 magnitude, we used these fields as 
zeropoint reference fields. 



We Stacked the post-Elixir processed, astrometrically cali- 
brated images using Swarp software by E. Berlin at the TER- 
APIX (Traitment Elementaire, Reduction et Analyse des Pix- 
els de megacam)'". By running Swarp, we created the mosaic 
image and the corresponding weight image (coverage map) 
whose pixel value is proportional to the exposure time. The 
reduced images were resampled, background subtracted, flux 
scaled, and finally coadded to one image. The interpolation 

http://terapix.iap.fr 



function we used in the resampUng isLANCZ0S3, which uses 
a moderately large kernel. The pixel scale of the resampled 
image is kept to the original MegaCam pixel size of 0"185. 
The astrometry transformation is handled by the Swarp pro- 
gram automatically, although the TNX WCS as derived in the 
previous section is not a FITS standard. 

For the background subtraction, we used the background 
mesh size of 128 x 128 pixels, and applied 3x3 filter box. 
The background mesh size and the bin size are chosen to bal- 
ance out the effect of bright stars (in general smaller than 128 



7 



1.0 



0.5 



0.0 



-0.5 



0.5 



^ 0.0 



-0.5 



l;8 



0.5 



0.0 



-0.5 



0.5 



0.0 



-0.5 



-1.0 



FLSll-cenlral 



FLS21-oenlral 



FLS32-FLS23 



FLS32-central 



■16 




FLS12-central 



FLS22- central 




FLS31 -central 



,0 19.0 20.0 21.0 22.0 23.0 24.0 
u* magnitude(AB) 



18.0 19.0 20.0 21.0 22.0 23.0 24.0 
u* magnitude(AB) 



Fig. 5. — After the zeropoint correction, the magnitude differences between fields are in the order of 0.01 magnitude. 



X 128 pixels) and the effect of the overall background gra- 
dient. We tried many background combinations to obtain the 
background as flat as possible and settled on the above param- 
eter values. The background-subtracted images are combined 
by taking the weighted average of each frame. The weight 
images are constructed by taking the bad pixel masks, and re- 
placing the bad pixel values to be 0, and the good pixel values 
to be 1 . The weight images constructed this way do not ac- 
count for small variation in the pixel response over the CCD 
(flat-field). Finally, the images are coadded using the follow- 
ing equation: 



In this equation, F is the value of a pixel in final coadded im- 
age, Wi is the weight value for the pixel, and fi is the value of 
the pixel in each individual image. The factor represents 
the flux scale for each individual image. The images are cal- 
ibrated to have the same photometric zeropoints after Elixir 
pipehne, but due to the airmass differences, the flux scales of 
the images are slightly different. In order to take into account 
the difference in zeropoints, we put the flux scale of the image 
having the least extinction to be 1.0. The flux scales of other 
images were calculated according to the equation: 

where k is the coefficient for airmass term and Xq and X, are 
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airmass values of the reference image and the image in con- 
sideration, respectively. The coadded gain of the final stacked 
image varies with position and is proportional to the number 
of frames used for producing each coadded pixel value. 

3.5. Photometric Calibration 

We performed photometric calibration on the mosaic im- 
age of each field using the photometric solution provided by 
the Elixir pipeline. During the production of mosaic images 
(§3.4), we rescaled each dither image so that each coadded 
image has the photometric zeropoint of the image with the 
least airmass value. Then, these magnitudes were corrected 
for the galactic extinction estimating the amount of galactic 
extinction from the extinction map of Schlegel et al (1998). 
Since the FLS field lies at moderately high galactic latitude, 
the amount of galactic extinction is relatively small. They are 
0.15 mag and 0.1 mag on average for the u* and ^-band re- 
spectively. In the next step, we derived an additional zeropoint 
correction necessary to account for some of the data that were 
taken under non-photometric condition (e.g., thin cirrus). 

To do so in M*-band, we used fields whose mosaic image is 
made of photometric dither images only. The fields, FLS21, 
FLS22, and FLS31 are such fields, and we used them as ref- 
erence photometry fields to derive photometric zeropoint cor- 
rection of other M*-band fields. Figure 4 shows the compar- 
ison of M*-band magnitudes (before correcting for the non- 
photometric data) between different fields that overlap with 
each other. For the comparison of magnitudes, we used the 
total magnitude (MAG_AUTO) from SExtractor (See Section 
§4). In Figure 4, we can see that M*-band magnitudes of ob- 
jects in FLS21, FLS22, and FLS31 fields are slightly brighter 
than the same objects in the central field. Also, objects in 
fields other than the reference photometry fields have u*- 
band magnitudes similar or fainter than those of the central 
field. The M*-magnitude offsets between the central field and 
FLS21/22/31 are 0.03 magnitude and the rms in the offset 
values is of order of < 0.01 magnitude. This confirms that 
the u-band magnitudes of FLS2 1/22/31 fields are the most 
reliable. The final M*-band zeropoint of each field are cal- 
culated using this strategy, and we present the zeropoints in 
Table 2. The comparison of M*-band magnitudes of overlap- 
ping objects after the photometric caUbration are presented 
in Figure 5. Based on the photometric consistency between 
FLS21/22/31, we estimate the accuracy of the M*-band photo- 
metric zeropoint to be of order of 0.01 magnitude. 

For g*-band mosaic, we preformed photometry on dither 
images that are taken under photometric condition, and de- 
rived a necessary zeropoint correction by comparing the pho- 
tometry of the mosaic image and the photometric reference 
image. The ^-band zeropoint of the central field derived this 
way is also given in Table 2. 

4. CATALOGS 

4. 1 . Object Detection and Photometry 

In the central 1 deg^ region (the 03B run), the stacked u*- 
and ^-band images were registered to the ^-band image. Both 
bands have the same pixel scale, so they can be easily regis- 
tered with the xregister task in IRAF. After the registra- 
tion, source catalogs were created using dual mode photome- 
try with SExtractor (Bertin & Amouts 1996). Specifically, ob- 
jects are detected in the ^-band image, and photometry is per- 
formed on both the g- and M*-band images based on these po- 
sitions. The weight image (coverage map) produced by Swarp 




Fig. 6. — The distribution of the values of CLASS_STAR parameter as 
a function of magnitude are presented. At magnitudes below flux limit, the 
stellarity index of the objects are distributed between 0.35 and 0.7. The fea- 
ture of the CLASS_STAR distribution at faint objects (< 27 mag) is due to 
the difficulties when calculating steUarity index with isophotal areas of the 
objects. 



was used as the WEIGHT_IMAGE. The LOCAL background is 
estimated using a 128 x 128 pixel background mesh with 3 
X 3 median boxcar filtering, identical to the default Swarp 
settings. 

We filter our images with a gaussian convolution kernel 
(gauss_4 . 0_7 x 7 . conv) matched to our average seeing 
conditions (FWHM^0."85). Finally, source detection was 
performed on the background-subtracted, filtered image by 
looking for pixel groups above the detection threshold. After 
playing with various parameter combinations and inspecting 
the results by eye, we settled on a minimum of 8 connected 
pixels and a 1 .2a minimum detection threshold. We also tried 
various values for the deblending parameter to separate ob- 
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TABLE 3 

Selected 15 entries of M*/g band source catalog 



obj ID 


RA(J2000) 


DEC(J2000) 


"auto 


€ 


11* 

"aper 


e 


Sauto 


e 


&WER 


e 






extg 


1001 


17:14:13.301 


59:15:19.199 


24.992 


0.051 


25.016 


0.042 


24.726 


0.066 


24.699 


0.052 


0.010 


0.120 


0.088 


1002 


17:17:11.895 


59:15:27.342 


24.882 


0.043 


25.200 


0.052 


24.834 


0.063 


25.112 


0.073 


0.006 


0.169 


0.125 


1003 


17:15:52.126 


59:15:26.820 


24.872 


0.040 


25.184 


0.048 


24.796 


0.053 


25.135 


0.064 


0.086 


0.151 


0.111 


1004 


17:13:19.966 


cn. 1C.1/C /111 

59:15:16.411 


24.955 


0.050 


25.408 


0.061 


24.582 


0.069 


25.144 


0.092 


0.024 


0.114 


0.084 


1005 


17:19:55.115 


59:15:21.739 


25.356 


0.080 


25.615 


0.070 


25.364 


0.064 


25.744 


0.063 


0.928 


0.141 


0.104 


1006 


17:14:53.745 


59:15:25.502 


25.628 


0.088 


26.223 


0.156 


25.666 


0.071 


26.027 


0.102 


0.586 


0.124 


0.092 


1007 


17:13:55.459 


59:15:21.272 


25.644 


0.091 


26.233 


0.123 


25.627 


0.072 


26.334 


0.109 


0.232 


0.118 


0.087 


1008 


17:18:35.120 


59:15:27.905 


26.073 


0.138 


27.508 


0.390 


26.164 


0.110 


28.560 


0.752 


0.725 


0.158 


0.116 


1009 


17:15:37.141 


59:15:28.592 


26.565 


0.189 


26.240 


0.122 


26.578 


0.181 


26.348 


0.128 


0.693 


0.146 


0.107 


1010 


17:15:33.977 


59:15:19.844 


23.060 


0.009 


24.161 


0.019 


22.738 


0.013 


23.741 


0.026 


0.028 


0.145 


0.106 


1011 


17:14:33.853 


59:15:19.858 


23.000 


0.008 


23.381 


0.010 


22.917 


0.011 


23.295 


0.013 


0.048 


0.121 


0.089 


1012 


17:17:52.075 


59:15:25.035 


23.869 


0.018 


24.589 


0.030 


23.858 


0.037 


24.674 


0.068 


0.029 


0.180 


0.132 


1013 


17:18:54.719 


59:15:23.401 


23.485 


0.013 


23.992 


0.016 


23.444 


0.015 


23.960 


0.019 


0.959 


0.148 


0.109 


1014 


17:16:11.917 


59:15:26.175 


23.960 


0.018 


24.495 


0.025 


23.890 


0.023 


24.382 


0.031 


0.036 


0.153 


0.113 


1015 


17:17:19.797 


59:15:29.375 


26.003 


0.113 


26.299 


0.139 


26.031 


0.169 


26.255 


0.194 


0.131 


0.172 


0.127 



^The AUTO magnitudes are calculated by SExtractor, using BCron-like elliptical aperture (Kron 1980). 
''The APER magnitudes are calculated over a circle with 3" diameter 

■^The stellarity represents CLASS_STAR calculated with SExtractor. The values are distributed between (galaxy) and 1 (star). 

*^The galactic extinction values are calculated using the extinction map of Schlegel et al.(1998) 

jects that are close together, and we adopted a deblending 
threshold of 32 and a deblending minimal contrast of 0.005. 

Using dual mode photometry, we detect ^^200,000 sources 
in the central 1 deg^ field. We also performed a single mode 
photometry on the M*-band image with the same configuration 
file, to include w-band objects that are not detected in g-band. 
About -^150 objects are detected in M*-band without ^-band 
detection. 

We applied a nearly identical method and parameters for 
source detection on the seven 04A FLS M*-band images. 
Since there are several very bright stars in the 04A M*-band 
images, we chose a bigger background mesh size - 256 x 
256 pixels to avoid the over- subtraction of the background 
near bright stars. Since the seeing during 04A was slightly 
worse than that of 03B (Table 2), we adopted a slightly larger 
gaussian convolution kernel (gauss_5 . 0_9 x 9 . conv). As 
stated earlier, the gain for each stacked image is determined 
to be the original gain times the number of stacked images. 

Due to interchip gaps and our 30~40" dither steps, the ef- 
fective exposure time varies from pixel-to-pixel. We there- 
fore applied the Swarp weight image when performing pho- 
tometry. Through the visual inspection of the weight image, 
we determined the value of WEIGHT_THRESH parameter for 
weighting in SExtractor. The value of the pixels in interchip 
gaps are measured, and used as threshold parameter value. 
This procedure prevents the rise of unusual sequence in the 
magnitude vs. magnitude error plot. 

In the catalogs, we present the apparent magnitude of the 
objects in two different ways (aperture & total magnitude). 
We measured the aperture magnitude using 3" diameter aper- 
tures. We adopted the auto magnitude (MAG_AUTO) from 
SExtractor as the total magnitude, which is calculated using 
the Kron elliptical aperture (Kron 1980). The related param- 
eters, PHOT_AUTOP ARAMS were detemiined to be 2.5, and 
3.5, which are the values in default configuration file. 

4.2. Star-Galaxy Separation 

We use the SExtractor stellarity index to separate star-like 
objects from extended sources. SExtractor uses 8 measured 
isophotal areas, peak intensity and seeing information to cal- 
culate the stellarity index. With the SEEING_FWHM parame- 
ter value and the neural network file as the inputs, SExtractor 
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Fig. 7. — Comparison between the positions of matched objects in our 
g-band catalogs and SDSS g-band catalogs. Offsets are computed as SDSS 
minus CFHT positions. 



gives the stellarity index (CLASS_STAR), which has values 
between 1 (starlike object) and as the result of the classifi- 
cation. Figure 6 shows the distribution of stellarity index as 
a function of magnitude in the central 1 deg^ field. We ran- 
domly picked objects with different apparent magnitudes and 
stellarity index values, and visually inspected them to see how 
useful the stellartiy index is for the point source vs. extended 
source classification. Our visual inspection reveals that bright 
objects (^ < 21 mag, and u* < 20 mag) with the stellarity in- 
dex > 0.8 are found to be stars. Through the visual inspection, 
we adopted the following criteria for separating point sources 
from extended objects: 

CLASS_STAR>0.8 for 20< U <23 
CLASS_STAR>0.8 for 21< g <23 
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Fig. 8. — Comparison between the positions of matclied objects in our w*-band catalogs and i?-band, IRAC catalogs. The offsets are computed as other survey 
minus CFHT. 



In addition to these criteria, we weed out brigiit saturated 
stars (g <21 mag, u* < 20 mag), based on SDSS point source 
catalogs. The subtraction of objects matched with SDSS stars 
is a good method for removing stars from our catalogs. 

On the other hand, it is difficult to address the nature of faint 
(g > 23 mag, u* > 23 mag) objects with the stellarity value 
alone, because faint, distant small galaxies are hardly resolved 
under 1" seeing conditions. Fortunately, number counts tend 
to be dominated by galaxies at m*,^ > 23 mag. Therefore we 
do not attempt to separate stars from galaxies at these flux 
levels. Also, the stellarity index distribution shows a converg- 
ing feature between 0.35 and 0.7 at magnitudes below flux 
limit (g > 26.5 mag, u* > 26.2 mag). This is thought to be 
the effect of difficulties in determining isophotal areas of faint 
objects. 

4.3. Catalog Formats 

As an example, 15 entries of the central 1 deg^ u*- and ^- 
band merged catalog is presented in Table 3. All magnitudes 
are given in the AB magnitude system. We include the to- 
tal magnitude (MAG_AUTO) and the aperture magnitude (3" 
diameter) in the catalog. The magnitudes are not corrected 
for the galactic extinction, but the extinction values are pre- 
sented as a separate column. Magnitude errors are the outputs 
from SExtractor. Considering the error in zeropoint and the 
extinction correction, we believe that the minimal error in the 
magnitude will be about 0.02 mag. Stellarity indices are cal- 
culated in the ^-band image for the central field. 

5. PROPERTIES OF THE DATA 

5.1. Astrometry 

We derived the final astrometric solution using SDSS 
sources (both point and extended sources). The rms error be- 
tween SDSS point sources and our sources are calculated to 
be Q'.'l<dr <0f'15. Considering the absolute astrometric er- 
ror of 0"1 in SDSS (Pier et al. 2003), our catalogs are fliought 
to have the positional error of roughly 0."15 - 0."2. Figure 7 
shows the positional differences between SDSS sources and 
our sources over the entire mosaic of the central 1 deg^ field. 



To address the accuracy of astrometric calibration in other 
ways, we compared our source positions with other avail- 
able catalogs for the Spitzer FLS. The catalogs compared are 
KPNO R-bmd data (Fadda et al. 2004) and flie Spitzer IRAC 
data (Lacy et al. 2005). The results are given in Figure 8. The 
average rms error, dr = \J {lS.acosS^ + IS.S^'), is about '~0."28 
with respect to /?-band data. /?-band positions have a sys- 
tematic offset of Aa = 0f'09 ± 0."27 and = 0."19 ± 0^'31 
with respect to CFHT positions. The systematic offset is orig- 
inated from the fact that the reference catalog for 7?-band data 
is GSC II stars, which is known to have a systematic offset 
with respect to SDSS reference positions (Fadda et al. 2004). 
Therefore, we consider the astrometry difference between R- 
band and CFHT data is well understood. For IRAC sources, 
the average rms error is dr ~ 0"25. The offset of IRAC 
positions according to CFHT positions are calculated to be 
Aa = -0."08 ± 0f'40 and AJ = 0."04 ± Of'38. The FLS IRAC 
data have positional accuracy of 0"25 rms with respect to 
2MASS (Lacy et al. 2005). The difference between 2MASS 
positions and SDSS positions caused the offset between IRAC 
positions and CFHT positions. The results add support to the 
astrometric accuracy of our dataset. 

5.2. Photometry Transformation 

As we have shown in Figure 1 , MegaCam filters are slightly 
redder than SDSS counterparts. Therefore, we examined 
the correlations between u* magnitudes and SDSS u magni- 
tudes. To do so, we matched non-saturated point sources with 
18.5< u <20 from our catalog with the SDSS point source 
catalog, and compared the differences in their magnitudes. 
Figure 9 shows the comparison, and it demonstrates that the 
difference between the m* and u varies from 0. 1 to 0.6 magni- 
tude as a function of the color of the object. We fitted the rela- 
tion with a first order polynomial to construct the conversion 
formula between CFHT u* and SDSS M-band magnitudes, and 
obtained the following relation from the linear least-square 
fitting. Although there exist some scatter, the conversion for- 
mula could be a good reference in use of CFHT u* magnitude. 

^mag{USDSS-UCFHT) = 
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Fig. 9. — Empirical difference between CFHT «*-band and SDSS u-band 
magnitudes as a function of {u-g) color. We used least-square fit to find the 
conversion equation from one M-band to another. 



Fig. 10. — Calculated u* -u for various spectral types of stars. The stellar 
templates from 05V to F9V stars were used in calculation, since in case of 
red stars, M-band is not good filter to sample their red colors. 



(-0.081 ±0.006) + (0.237 ± 0.004) x (u-g)sDss 

We checked the above empirical relation between Mega- 
Cam u* and SDSS u by calculating (u*-u) vs. (u-g) using 
the filter curves in Figure 1 and the stellar templates from 05 V 
to F9V stars (Silva & Cornell 1992), or the empirical galaxy 
templates at different redshift (Coleman et al. 1980). The the- 
oretically calculated relations of {u*-u) vs. (u-g) agree well 
with the empirical relations (Figure 10, 11). 

On the other hand, the difference between SDSS g-band 
magnitude and CFHT g-band magnitude is considerably 
small. Figure 12 shows their difference along (u-g) color. 
Compared to Figure 9, the magnitude difference is small (they 
scatters within 0.15 mag), but the differences still have a ten- 
dency of increasing toward red color. We applied least-square 
fit to g-band magnitude difference to find the following con- 
version relation. 



Armg(gsDSS-8cFHT) = 
(-0.061 ±0.007) + (0.057 ± 0.004) x (u-g)sDSS 

5.3. Galaxy Number Counts 

To show the depth and homogeneity of the survey, we 
present the galaxy number counts in m* and g-band. The num- 
ber counts are constructed for both the deep central 1 deg^ 
field as well as the shallower M*-band coverage of the entire 
FLS fields. Then, the results are compared with galaxy counts 
from other studies. The direct comparison is possible for the 
counts constructed from SDSS filters (Yasuda et al. 2001), 
while, for other studies using lohnson U, we converted the 
Johnson {/-band magnitude to our u* magnitude using the 
M-f / ~ 0.8 conversion relation of galaxy colors from Fukugita 
et al. (1995), and then applying the mean of our empirical 
conversion relation between u* and u(u*-u = 0.4 mag). For 
g magnitude, we used 5-band number counts from the litera- 
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Fig. 11. — Calculated difference between the CFHT «*-band and SDSS 
M-band magnitudes as a function of (w - g) color. The galaxy templates used 
are elliptical, spiral, & irregulars from Coleman, Wu. & Weedman 1980. We 
redshifted the galaxies to various redshifts, and calculated their colors using 
filter curves. 

tures (Metcalfe et al. 2000; Capak et al. 2004; Amouts et al. 
2001; Huang et al. 2001), converting mag to g mag accord- 
ing to the relation of g-B = 0.3 ~ 0.6 (Fukugita et al. 1995). 
The conversion relation from Johnson U to u, and from Btog 
filter varies according to the shape of the galaxy spectral en- 
ergy distribution, so we adopt the mean difference to translate 
one filter to another. 

The number counts are shown in Figure 13. When con- 
structing the number counts, we excluded stellar objects us- 
ing the stellarity cut introduced in section §4.2. For objects 
brighter than g <2l mag and u* < 20 mag, we matched the 
SDSS point source catalogs with our catalogs, and subtracted 
matched objects which are considered to be stars. Error bars 
on the figure represent the poisson errors only. Also note that 
our counts are not corrected for completeness. 

With the exception of SDSS studies (Yasuda et al. 2001), 
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Fig. 12. — Empirical difference of dfTT g-band magnitude and SDSS 
g-band magnitude according to (u-g) color. Because CITiT filter is redder 
than SDSS filter, there is also a small tendency of magnitude difference along 
the color. But it is not difficult to say that the difference is negfigible. 

galaxy counts are typically derived from deep, but small area 
surveys (~0.2-0.3 deg^). Our image covers from 0.94 deg^ 
(in case of g-band) to 5.13 deg^ (in case of shallower m*- 
band). For u* < 24 mag, we include the shallower M*-band 
data to maximize the area coverage. We can say that our 
number counts are better than other number counts in terms 
of the area coverage. For future use, we present a table that 
summarizes the result of galaxy number counts (Table 4). To 
construct the number counts, we use the auto magnitude in 
SExtmctor. Our number counts are in good agreement with 
the counts from other studies to u* ~ 24.8 mag, ^ ~ 25.2 mag. 
Beyond u* ~ 24.8 mag and g ^ 25.2 mag, our data starts to 
tail off, and this shows that our u* and g catalogs are nearly 
100% complete down to the above magnitude limits. These 
limiting magnitudes have the uncertainty of 0.3 magnitudes 
due to the fluctuation in the galaxy number counts at the faint 
end. We also obtained the number counts for the shallower 
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Fig. 13. — Galaxy number counts in FLS region from our u", g-band 
images. Results from other studies are overplotted. Magnitudes are in AB 
system. 

M*-band data, the results are illustrated in Figure 14. In case 
of the deepest field among the shallower M*-band, the 100% 
limit is 24.2 magnitude. 

5.3.1. Completeness 

To inspect the completeness of our survey in more detail, 
we used two independent methods. First, we compared our 
galaxy counts with the deeper galaxy counts from other stud- 
ies. Second, we made artificial images having the same prop- 
erties with our observation (seeing, crowding, gain & back- 
ground) using art data package in IRAF. The artificial im- 
age contained both point sources and extended sources. The 
extended sources are defined to have 40 % of whole galax- 
ies as elUptical galaxies and the remainders for disk galaxies, 
with minimum redshift of z ~ 0.01. Poisson noise is added 
to the background. Then, we performed photometry with 
the same configuration file we used to make catalogs of ob- 
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Fig. 14. — Galaxy number counts in FLS region from the shallower u*- 
band images. Different solid lines indicate the different fields. The shallow- 
est field, [FLS23] shows 100% completeness to about 23.6 magnitude. The 
deepest field, [FLSll] is 100% complete down to 24.2 mag. The dashed line 
is the result from the central 1 deg^ region. 
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Fig. 15. — The completeness test. Completeness was inspected by dividing 
galaxy number counts with median number counts from other surveys (solid 
line). The dashed fine shows the estimation of completeness from artificial 
image made. 

served sources. The results are shown in Figure 15, where 
the solid fine and points indicate the completeness from the 
first method, and the dashed Une represents the completeness 

estimated from the simulation. Figure 15 shows that the data 
is 100% complete down to ^ ~ 25.2 mag and u ~ 24.8 mag, 
agreeing the value found in the previous section, and is ^50% 
complete to 26.5 magnitude for g-band, 26.2 magnitude for 
M*-band. 

For the shallower M*-band data, we examined the complete- 
ness using the first method (See Figure 16). The 50% com- 
plete limiting magnitudes are 24.8-25.4 mag, which are con- 



14 





' 1 ' 








1 < < < 1 , , , 1 < < 







U"* magnitude(AB) 



Fig. 16. — Completeness of the survey as a function of magnitude for the 
shallower M*-band. The dashed line is the completeness of the central 1 deg^ 
region M*-band. 

sistent with the estimated value in Table 2. 

5.3.2. Reliability 

We also examined the reliability of our detections. We cre- 
ated a negative image by multiplying -1 to the mosaic image 
and performed photometry on the negative image. The num- 
ber of objects detected in the negative image gives us an esti- 
mate of how many false detections exist in our catalogs. The 
false detection rate is defined as 

Nfaise(m,m+Am) 

N(m,m+ Am) 

where Nfaisei'n,m+ Am) is the number of objects detected in 
the negative image over the magnitude bin of (m, m+Am), and 
N{m,m+Am) is the number of objects detected in the original 
mosaic over the same magnitude bin. The false detection rate 
is 0.5, when there is only noise in the image. 

The number of spurious detections are relatively small (less 
than 0.7%) until g ~ 26.5 mag and u* ~ 26.0 mag, increases 
steeply toward fainter magnitudes. The presence of false de- 
tections at brighter magnitudes (which is above the Umiting 
magnitude) is due to the effects of streaks around bright stars. 
On the other hand, the majority of false detections at fainter 
magnitudes is generated from the noise. The similar tendency 
is emerging on the shallower M*-band images. The false de- 
tection rate remains under 1.8% until the Umiting magnitudes 
for the shallower M*-band. 
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6. SUMMARY 

The entire ~5 deg^ region of the Spitzer First Look Survey 
area have been observed with the MegaCam on the CFHT 3.6 
m telescope. From the data, we created nine final mosaic im- 
ages, photometrically and astrometrically calibrated. Two of 
them are the deeper u*, g-band images on the central 1 deg^, 
and the remaining seven are the shallower M*-band images 
over the whole 5 deg^. The average seeing for the central 1 
deg2 field is 0"85, while the seeing is ~1."00 otherwise. The 
central field goes as deep as g ~ 26.5 mag and u* 2± 26.2 mag 
at 5a flux hmit within 3" aperture (or ~50% completeness), 
while the shallower M*-band data goes to 25.4 mag. Our data 
are deep enough to detect optical counterparts of many IRAC 
and MIPS sources in the FLS region. We recalibrated the as- 
trometric solution, and the final absolute astrometric uncer- 
tainty is at the level of 0." 1 when compared with SDSS point 
sources. With other surveys on FLS, the rms error is under 
0"3, and it reflects merely the astrometric accuracy of the 
other datasets. The high astrometirc accuracy allows this data 
to be used efficiently for multislit or fiber spectroscopic ob- 
servation of faint galaxies in the FLS. Using dual mode pho- 
tometry, with g-hmd as the reference image, about 200,000 
extragalactic objects are detected in our central 1 deg^ (ef- 
fective observed area reaches ~ 0.94 deg^) down to g ~ 26.5 
mag. There are ^50,000 sources per deg^ in the shallower 
M*-band images down to u* ~ 25.4 mag. Our galaxy number 
counts are consistent with other studies, and the number of 
spurious detections is less than 0.7% below our Umiting mag- 
nitudes. Our data could be used in many ways: for example, 
i) selection of high redshift objects using broad-band dropout 
technique, ii) improvement of photometric redshifts, and iii) 
investigation on rest-frame UV properties of high redshift in- 
frared sources. These studies are currently under way. The 
mosaic images and catalogs will be available through a public 
website in near future. 
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TABLE 4 

The result of galaxy number counts 
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